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ABSTRACT 

We have analysed the 24-yum properties of a radio-selected sample in the Subaru- 
XMM/Newton Deep Field in order to explore the behaviour of the FIR/radio relation at high 
redshifts. Statistically, the correlation is described by g24, the ratio between the observed flux 
densities at 24 /jm and 1 .4 GHz, respectively. Using 24-/im data results in considerably more 
scatter in the correlation than previous work using data at 60-70/im. Nevertheless, we do 
observe a steady correlation as a function of redshift, up to z ~ 3.5, suggesting its validity 
back to primeval times. We find #24 = 0.30 + 0.56 for the observed and ^24 = 0.71 + 0.47 
for the ^-corrected radio sample, based on sources with 300 ^/Jy < S1.4GHZ < 3.2 mJy and 
24yum detections. A suitable ^-correction given by a M82-like mid-IR template suggests no 
extreme silicate absorption in the bulk of our radio sample. Using thresholds in ^24 to identify 
radio-excess sources, we have been able to characterise the transition from radio-loud AGN to 
star-forming galaxies and radio-quiet AGN at faint (< 1 mJy) radio flux densities. Our results 
are in broad agreement with previous studies which show a dominant radio-loud AGN popu- 
lation at > 1 mJy. The rest-frame U - B colours of the expected radio-excess population have 
redder distribution than those that follow the correlation. This is therefore a promising way to 
select obscured Type-2 AGN, with a radio loud nature, missed by deep X-ray observations. 
Spectroscopic follow-up of these sources is required to fully test this method. 

Key words: starburst - galaxies: active - galaxies: high-redshift - galaxies: classification, 
colours - galaxies. 



1 INTRODUCTION 

The far-infrared/radio (FIR/radio) relation is one of the tightest cor- 
relations known in astronomy. Established over three decades ago 
(e.g. Van der Kruil 1971; Helou et al. 1985), it covers about five or- 
ders of magnitude in luminosity (Condon lll992l; Garret fc002h . The 
correlation is tightest for local star-forming galaxies and relates the 
integrated thermal emission, from the dust present in star-forming 
regions, with the non-thermal synchrotron emission produced by 
relativistic particles (cosmic rays) accelerated by supernovae explo- 
sions (Harwit & Pacini Tl975l) . A complete explanation of the origin 
of this remarkable correlation is essential for our understanding of 
how galaxies formed and evolved, particularly if the correlation ex- 
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tends from primeval galaxies to the fully mature g alaxies observed 
locally, as recent studies have hinted (Kovacs et al 2006). 

Being relatively direct probes of recent star formation, FIR 
and radio data are amongst the most powerful tools available to test 
whether primordial galaxies shared the same properties as galaxies 
present today. Low-frequency radio data has been used to estimate 
red shifts using radio-to-submm flux density ratios (e.g. Carilli & 
Yuri 2000), and to determin e the co smic star formation history (e.g. 
Haarsm dboool ; Ivison et al ll2007ah . facilitating the investigation of 
the cosmic IR background (CIRB). If the relation between flux den- 
sities observed in the FIR and radio wavebands holds at high red- 
shifts, then we can exploit ultra-deep high-resolution radio imag- 
ing to resolve the CIRB and to explore the properties of luminous 
distant galaxies - fo r exam ple, studies of subm m (or "SCUBA") 
galaxies (Smail et al j|l997l ; Hughes et al ]|l998l) which are encum- 
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bered by the confusion encountered in the submm (rest-frame FIR) 
wave band, and studies of active galactic nuclei (AGN; Beelen et al. 
l2006h that are subject to selection biases caused by obscuration at 
shorter wavelengths. 

While the combination of Spitzer observations at 24 and 70 pm 
and radio observations at 1.4 GHz support the universality of the 
FIR/radio correlation up to a redshift of approximately unity (Ap- 
pleton et al. ll2004) . extending this correlation to higher redshifts 
requires extremely deep and complementary observations at radio 
and IR wavelengths, as well as reliable estimates of redshift. The 
work presented here is based on a deep multi-wavelength study in 
the Subaiu-XMM/Newton Deep Field (SXDF). We have exploited 
1.4-GHz images from the Very Large Array (VLA), reaching a 5<r 
flux density threshold of ~ 35 /jJy and covering an area of approx- 
imately 1 deg 2 . The region has also been observed to impressive 
depths at optical (SuprimeCAM; Furu sawa e t al. in preparation), 
near-IR (UKIDSS-UDS; Lawrence et al ll2007l) and mid-IR ( Spitzer 
Wide-area IR Extragalactic Survey - SWIRE; Lonsdale et al l2003l) 
wavelengths. 

Starburst galaxies, radio-quiet and radio-loud AGN are the 
three main populations which can be observed to high redshifts 
with the data employed in this work. It is known that radio-quiet 
AGN (i.e. Seyfert galaxies) follow the FIR/radio correlation (Roy 
et al ll998l) . but radio-loud AGN may vary considerably in radio lu- 
minosities due to their extra core and jet radio emission, thus weak- 
ening and biasing the correlation. In theory, it is possible to identify 
these radio-loud AGN via their deviation from the expected relation 
and therefore to determine the fraction of these so-called "radio- 
excess" sources as a function of radio flux density. This would then 
allow us to constrain the predictions of synthesis population mod- 
els for the transition at sub-mJy radio fluxes from radio-l oud AGN 
to starbu rst/rad io-quiet AGN (Dunlop & Peacock lll990l : Jarvis & 
Rawling jliooj) . 

The aim of this work is to explore and discuss the correlation 
between the monochromatic fluxes at 24yum and 1.4GHz across a 
broad range of redshifts and test the capabilities of the relation for 
selecting radio-loud AGN (e.g.. Donley et al 112004!) . 

In Section [2] we present the multi- wavelength data available 
for the SXDF field and used in this work. Section|3]shows the ca- 
pabilities of the data to extend the 24yum/1.4GHz correlation to 
higher redshifts and Section [4] describes the process of detecting 
radio-loud AGN at sub-mJy 1.4-GHz flux densities. The discus- 
sion and conclusions are presented in Section[5]and[6]respectively. 
Throughout this paper we have us ed £l m = 0.3, Aa = 0.7 and 
H = 70 km s -1 Mpc" 1 (Spergel et al ll2007l) . 



2 SXDF OBSERVATIONS 

2.1 Radio observations at 1.4 GHz 

Radio o bservations of the SXDF were carried out by Ivison et al. 
M2007bl) during July 2003 using the VLA at 1.4 GHz. The total 
integration time was approximately 60 hours, reaching an r.m.s. 
depth of 7 /jJy beam -1 near the centre of the field. The radio map is 
centred at RA. 2h 18m, Dec. -5° 0', with a synthesised beam of 
1.86arcsec x 1.61 arcsec at position angle 15°. The source detec- 
tion was carried out using a 5<x threshold, sufficient to avoid signif- 
icant contamination by spurious sources, obtaining a total number 
of 563 radio sources in an area of ~0.45 deg 2 . 

In order to exploit the larger area covered by previous multi- 
wavelength observations in the SXDF, we have also used the 1.4- 



GHz catalogue presented by Simpson et al. H20061) . Their obser- 
vations comprised thirteen VLA pointings reaching a mosaic r.m.s. 
noise of ~ 1 2-20 /iJy beam -1 . The Simpson et al. catalogue contains 
505 sources and covers 0.8 deg 2 to a flux density limit of lOOyuJy. 

By combining both radio samples, we obtained a total sam- 
ple of 828 radio sources. Extended sources and obvious doubles, 
identified by eye, were considered as single emitters. The distribu- 
tion of radio flux densities peaks at about 100//Jy and reaches flux 
densities as faint as 35 p]y (5tr) and as bright as 80 mjy. 



2.2 Spitzer observations 

The SXDF was observed by Spitzer as part of the Spitzer Wide - 
area Infrared Extragalactic Survey (SWIRE; Lonsdale et al 11200 3h 
Legacy programme. The SXDF is included in one of the seven 
high-latitude SWIRE fields, the XMM-LSS field, with coverage of 
9.1 deg 2 . Observations were obtained at 3.6, 4.5, 5.8 and 8.0pm 
with IRAC and at 24, 70 and 160 yum with MIPS (Surace et al. 
2005). In order to explore the FIR/radio correlation, we have used 
only the detections at 24 pm. It is known that use of the monochro- 
matic 24-pm flux de nsities yields a larger scatter than the 70-yum 
data (Appleton et al. 2004); however, we are limited by the rela- 
tively shallow coverage at 70 pm and the more severe confusion at 
that wavelength. 

To explore the FIR/radio relation as a function of redshift, we 
use the combined radio source catalogue described in §2. 1 to cross- 
match it with the conservative (lOcr) SWIRE source catalogue re- 
trieved from the Spitzer Centre Archive. We obtained 370 coun- 
terparts with fluxes S^Vd ^ 0.4 mjy using a 12.5-arcsec-diameter 
aperture for the photometry. 

We then extended the catalogue to fainter 24-pm flux densities 
by measuring the observed flux density on the 24-pm map, using 
the same diameter aperture centred at the radio positions. This al- 
lowed the inclusion of 125 source s with 24-pm flux densities of 
S24fim 200yuJy (~4<t; Shupe et al. l2007h . In this paper, for all the 
remaining radio sources we used upper limits at Si^m = 200 ply. 

2.3 Subaru and UKIRT Observations 

The SXDF is home to the very deep optical survey undertaken 
by Subaru/SuprimeCAM (Miyazaki et al. 2002). This comprises 
five overlapping pointings, providing broad-band photometry in the 
BVRi'z' filters to typical 5-cr (AB magnitude) depths of B - 27.5, 
V = 26.7, R = 21.0, i' = 26.8 and z' = 25.9, respectively (2- 
arcsec-diameter apertures) . The s eeing in the composite images is 
-0.8 arcsec (Sekiguki et al 1120051 Furusawa et al j|2008h . 

The central region of the SXDF is also being observed at near- 
IR wavelengths with the Wide-Field Camera (WFCAM) on the 3.8- 
m United Kingdom Infrared Telescope (UKIRT). As the deepest 
tier of the U KIRT Infrared Deep Sky Survey (UKIDSS, Lawrence 
et al j|2007l) , the Ultra Deep Survey (UDS) aims to cover 0.8 deg 2 
to A" A b = 25. In this work, we use the UKIDSS First Data Release 
(DR1; Warren et al.|]2007) where the 5-cr point source depths (in 
AB magnitudes) at J and K are 23.61 and 23.55, with seeing of 
0.86 and 0.76 arcsec, respectively. 

We have used the /f-band as the reference to cross-match with 
all other near-IR and optical images. Unfortunately, the UDS map 
misses part of the radio coverage of the SXDF field, limiting the 
radio sample to an area with 726 radio sources, out of which 639 
(88 per cent) have clear counterparts in the UDS K data within a 
search radius of 1.5 arcsec. 
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2.4 Photometric redshifts 

We obtained reliable photometric redshifts for 586 of the radio 
sources detected at near-IR wavebands. The estimation was carried 
out using all available detections from the B,V,R,i',z',J,K, 3.6- 
yum and 4.5-yum photometric bands and considering a conservative 
criteria: not including either those sources with poor x 2 fitting nor 
those contaminated by halos from nearby stars. 

The spectral energy distribution (SED) fitting procedure pro- 
vides a photometric estimation consistent within Az/(1 + z) = 
0.05 ± 0.04 for the UDS galaxies in general (see Cirasuolo et al. 
1120071 for details), and spectroscopic redshif ts available for 63 of 
the 586 radio sources (Yamada et al. ll2005L Simpson et al. ll2006l 
and the NASA Extragalactic Database) indicate that the photomet- 
ric redshifts of the radio sources are also accurate at this level. The 
precision of this estimation is good enough to allow an adequate 
^-correction to the data (see Figure [T). The resulting redshift dis- 
tribution of the radio sources peaks at about z ~ 0.7 and includes 
galaxies up to z ~ 3.5. 

Note that since the photometric redshift estimation depends on 
the K-band detection, the sample is basically mass limited. Dwarf 
galaxies with faint radio emission may be missed in the K-band, but 
massive ones are expected to be detected up to very high redshift. 
Heavily obscured galaxies, however, are not drastically affected. 



S 1.4GHz > 3.2 mjy (10~ 25 Jy; large squares in Figure [T), charac- 
terised to be AGN-dominated by Condon Hl992l) . have qnn < -1, 
well below the bulk of the data points and similar to the estima- 
tion obtained from a composite radio-loud AGN template (Elvis 
et al II 19941) . Since we expect and want a correlation valid for star- 
forming systems, sources with S 14GHz > 3.2 mJy are excluded from 
our analysis of the FIR/radio correlation, basically because these 
sources are radio-loud AGN that introduce biases in the statistics. 

Using the biweight estimator^ (Beers et al. 1990) to charac- 
terise the data, we find a central location (mean) and scale param- 
eter (sigma; all ^24 values quoted in this paper will similarly be 
given in this way) given by q 2 4 = 0.66 ± 0.39 (based on all de- 
tections at 24yum including sources without estimated redshifts). 
The sub-sample of sources with < z < 1 gives an observed 
q 2 4 = 0. 80 + 0.33, in agreement with the results of Appleton et 
al. H2004I) . i.e. q M = 0.84 ± 0.28. The data also suggest a decreas- 
ing trend of q 2 4 as a function of redshift. A simple linear regression 
to the 24yum-detected sources with estimated redshift gives a de- 
pendency q 24 (z) = (0.85 ± 0.01) + (-0.20 ± 0.01) z. Nevertheless, 
this trend is probably seen because data have not been ^-corrected 
(Hogg et al ll2002l) as we describe later. It is worth noticing that the 
available 49 sources with spectroscopic redshifts support the valid- 
ity of the photometric estimation at high redshift (large circles in 
Figure [T}. 



3 THE FIR/RADIO RELATION 

Based on the analysis presented in Section [2] 369 sources (out of 
potentially 726 due to the missed UDS area) have both a 24-/jm 
detection with flux density S 24pm > 200 /jJy and a reliable redshift 
(49 of which are spectroscopic). 

We present in Figure Q] the FIR/radio relation in terms of the 
basic observable parameter, q 24 - the ratio between the monochro- 
matic fluxes at 24/im and 1.4GHz, log(S 24^/ S 1 4GHz) - as a func- 
tion of redshift. Ideally, the correlation should be calculated using 
the bolometric FIR flux, however the paucity of long-wavelength 
data makes this quantity difficult to estimate prior to the future ar- 
rival of SCUBA-2 (Holland et al. 2006) and/or the Herschel Space 
Observatory (Poglitsch et al. 2006). 



3.1 Non fc-corrected data 

Figure Q] shows that the vast majority of the sources detected at 
24 yum follow a tight correlation between ^24 and redshift, up to z ~ 
3.5. Upper limits have been included using 524,,,,, = 200 /jly and 
show probable evidence for a large number of radio-excess sources 
(i.e. radio-loud AGN) at the observed radio flux density regime. 

AT-corrections to the data have been tested assuming standard 
starburst galaxy templates for the 24^m detections and a simple 
power slope for the radio density fluxes. The expected variations 
of ^24 as a function of redshift are overplotted in Figure Q] and 
are used to study and characterise the distribution for star-forming 
galaxies. Nevertheless, the presence of radio-loud AGN sources in 
the sample introduces outliers in the correlation that may cover 
a wide range of ^24 values due to radio excess emission coming 
from core and jet structures. For example, the radio-loud AGN such 
as Cygnus A (a powerful double-lobed radio galaxy) and Centau- 
rus A (one of the closest known radio galaxies) have q 24 = -3.2 
and -1.1 values, respectively, far away from the main data trend 
from FigureQ] In particular, our data show that the population with 



3.2 fc-corrected data 

IR spectroscopic studies of local starburst galaxies have shown 
three main components in the mid-IR range (5-38/um): silicate 
bands around 10 and 18/jm, a large number of PAH emission fea- 
tures, and a slope of the spectral continuum, where an AGN may be 
present (Brandl et al. 2006). In order to model the ^-correction to 
the mid-IR emission from the bulk of the star-forming sources, we 
have convolved the SED of M82 (a local standard starburst galaxy 
with a distribution obtai ned fr om a fit to the observed photometry 
presented by Silva et al lll998l) . Arp 220 (obtained from photome- 
try listed in the NASA/IPAC Extragalactic Database and compiled 
by Pope et al.|[2006) and a model of an ultraluminous IR galaxy 
(ULIRG; L IR = 10 12 L Q ) given by Lagache et al. H2004I) . with the 
24-yum filter profile from SpitzerfMIPS. The radio emission is as- 
sumed to be well represented by a power law with a spectral index 
a = -0.7 , based on an average steep-spectrum radio source (Con- 
dor lll992l) . 

Note that assumptions of different mid-IR templates may 
largely change the ^24 values at high redshift, mainly due to the 
10-yum silicate feature redshifted to z ~ 1.5, where variations of 
even an order of magnitude in fc-correction are seen. On the other 
hand, if we consider an extreme error i n the ra dio spectral slope of 
a(a) = 0.5 (see findings in Garn et al j|2007l) . we find that the k- 
correction may change in ~ 0.3 dex at z ~ 3.5, a value much lower 
than the 24 /urn uncertainties. 

The overplotted lines, in Figure Q] show the expected change 
in the observed q 2 4 parameter as a function of redshift based on the 
different adopted templates. Apart from the composite radio-loud 
spectra, the ^-corrections have been normalised to the best fit to the 
data. 



1 In order to statistically describe the data, we used the biweight estimator 
which is resistant to outliers and is also robust for a broad range of non- 
Gaussian distributions. This is essential to remove remaining AGN from 
the statistics of the sample. 
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Figure 1. Black dots: ratio between the observed monochromatic fluxes at 1.4 GHz and 24 fjm, q 2 4 = log(S 2%im/S l.4GHzX as a function of redshift. Downward 
arrows: upper limits considering a mid-IR flux threshold given by 200 /iJy. Note that these data have not been ^-corrected. Sources without photo-z's are plotted 
at z = -0.1. Large circles: sources with spectroscopic redshifts. Large squares: sources with radio flux densities 5i.4ghz > 3.2mjy. Uncer tainties in the 
observed flux densities give an average ^24 mean error of ~ 0.1 dex, and error bars in redshift are approximately cr. a 0.1 (Cirasuolo at al 2007). Overplotted 
lines show different ir-correction factors, A'(24^m)/A'(1.4GHz), for the observed q 2 4 values as a function of redshift. The assumed mid-IR templates were 
convolved with the 24-/.im filter from Spitzer/MWS and are used to show the expected variation of the 524 values as a function of redshift. Three of these 
mid-IR templates are based on star-forming galaxies: two local standard massive starburst galaxies, M82 and Arp 220, and a SED model for an ULIRG galaxy 
with Lir = 10 L Q , as described by Lagache et al. 12004). All these three templates are assumed to have a simple spectral index a = -0.7 for th e radio 
emission, where S v K v", and normalised to the best fit to the data detections. Finally, a composite spectrum for radio-loud AGN given by Elvis et al. Ill994l) 
is also included, although as indicated by the locations of Cen-A and Cygnus-A, radio-loud AGN may be distributed over a wide range of q 2 4 values. In 
particular, this template is not normalised to the best fit but uses values from the actual composite spectrum at both wavelengths. 



A descending trend in q 24 as a function of redshift is expected 
from all three starburst-like ^-corrections, in agreement with the 
observed slope. The ^-corrected relation based on sources with 
mid-IR detections (> 200 yuJy) and estimated redshifts, changes to 
q 24 = 0.99±0.32, 1.45±0.38, and 1.08±0.37 (biweight estimator) 
using the M82-like, Arp 220-like and the Lagache et al. model, re- 
spectively. The M82-like fc-corrected q 24 values are presented in the 
upper panel of Figure[2]to demonstrate the improvement in the cor- 
relation as a function of redshift. Figure [3] shows a redshift-binned 
version of this plot (for the 24/mi detected sources), demonstrat- 
ing that both the mean and scatter of the M82-like ^-corrected 1724 
value remain constant out to z ~ 3.5; it also demonstrates that k- 
correcting with an Arp220-like spectrum leads to large variations 
in the relation with redshift. Based upon the behaviour of the M82- 
like ^-corrected distribution, and particularly the lack of a clear ex- 
cess in scatter for sources between 1 < z < 2, we suggest no strong 
MIR-silicate absorption for the bulk of the sources. 



The M82-like ^-correction to all sources with detections at 
24/jm in the range < z < 1 gives 24 = 1.03 ± 0.31, in agreement 
with previous Appleton et al. 120041) estimation (g 2 4 = 1.00±0.27). 
A simple linear regression to the M82-like ^-corrected data is given 
by 924(2) = (0.94±0.01)+(-0.01±0.01) z, suggesting a constant de- 
pendency (at lcr level) for the q 24 fraction with redshift. This result 
assumes that all galaxies in the sample can be fitted with a single 
M82-like template, although the fact that the correlation naturally 
appears to have the same mean and scatter as a function of redshift, 
up to z ~ 3.5, is highly suggestive that the FIR/radio correlation 
holds all the way back to primordial times. 

3.3 The correlation based on a complete radio sample at 
24yum 

The mean and scatter values quoted for the correlations described 
above may be being biased by incompleteness at 24 fim for sources 
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Figure 2. S"-corrected q 2 4 values as a function of redshift. The A-correction 
uses a M82-like template (Silva et al. 1998) for the mid-IR and a spectral 
slope a = -0.7 for the radio emission (S v K v°). Top: Full radio sample. 
Bottom: Radio sources with 1.4-GHz flux densities above 300 /jjy. Upper 
limits are considered at 5 24™ = 200 /jjy. 
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Figure 3. Binned mean and scatter of the M82-like ^-corrected 
24/jm/1.4GHz correlation (thick black data) as a function of redshift, based 
on the full radio sample. The overplotted lines correspond to the mean and 
one sigma limits, for the all-redshift distribution. Also shown (in thin black 
data) are the equivalent values for an Arp220-like ^-correction (shifted in 
redshift by an small factor in order to improve the visualisation only), which 
results in considerably more residual redshift variation. Note that the signif- 
icantly higher q 2 4 values for the Arp220-like ^-correction results from the 
upturn in the Arp220-like correction at z < 0.3 and due to the particularly 
extreme silicate absorption features (cf. Figure[T}. 




Swam MJy 

Figure 4. Mean q 2 4 value as a function of cut-off in radio flux density. 
Diamonds and crosses are based on the observed and M82-like ^-corrected 
data, respectively, both showing similar dependencies with the cut-off in ra- 
dio. The error bars represent the uncertainty on the mean value estimated us- 
ing boot-strap re-sample technique. The vertical dashed lines are at lOO^iJy, 
200 /iJy and 300 /iJy. Since sources with Si.4ghz > 3.2mjy (10~ 25 Jy) are 
essentially all radio-loud AGN (Condon 1 9841) . these were not included in 
the estimation for the mean value. 



with faint radio fluxes: a large number of star-forming dominated 
radio sources with S i.4ghz 5: 300 yuJy are not expected to have been 
detected at 24 jim. To demonstrate this problem, we plot in Figure|4] 
the mean ^24 as a function of a cut-off in the 1 .4-GHz flux density. It 
clearly shows a decreasing trend for the mean q 2 4 from fainter radio 
fluxes up to S 1.4GHz ~ 300/jJy, where a flattening starts to appear. 
This means that considering a high enough radio cut-off, the cor- 
relation becomes independent of the radio sample. We, therefore, 
base the correlation on a complete radio sample by constraining 
the number of sources to those with 5i 4GHz > 300 yujy. This radio 
sample is essentially complete at 24yum for star-forming galaxies, 
and all residual 24yum non-detection are associated with radio-loud 
AGN activity (see below). 

Using the biweight estimator, we find that for the radio pop- 
ulation with 5i 4GHz > 300yuJy and mid-IR detections (S^Afim > 
200 yuJy), the observed q 2 4 values has significantly lower mean and 
larger scatter for the (non-i-corrected) correlation q 2 4 = 0.30±0.56 
(it changes to ^24 = 0.70 ± 0.63 for the sub-sample of sources at 
0<z<l). 

A histogram of q 2 4 for the sources with 5i.4ghz > 300yuJy is 
shown in FigurefJ] It shows that the bulk of the undetected sources 
is found below q 2 4 = -0.18 due to the 200-/jJy threshold adopted 
for the mid-IR image. This value of 1724 implies sources three times 
brighter in radio luminosity than those from our observed mean 
value (qu) = 0.30, and ~ 10.5 times brighter than the observed 
me an val ue found by Appleton et al. 12004), <i?24> = 0.84. Yun et 
al. H200ll) have defined radio-excess sources as those whose radio 
luminosity is > 5 times the value predicted from the FIR/radio cor- 
relation. In this sense, radio-loud AGN based on 24-yum detections 
are sensitive to the assume d mean value for th e FIR/radio correla- 
tion (Appleton et al j|2004 Donley et al j|2005l Boyle et al j|2007t) . 
We find that 62 per cent of the radio sources (24 detections and 82 
upper limits) plotted in Figure [5] have 1^24 < -0.18. These values 
clearly suggest a large number of radio-loud AGN in this radio flux 
regime. 

It is worth noting that the incompleteness produced by 
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Table 1. Summary of the results presented in Section[3] The linear regressions are based on the detected sources at 1.4GHz and 24/im, using intrinsic error 
bars from the source extraction processes. 
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Figure 5. Histogram of the 924 values for the radio sources with radio 
flux densities brighter than 300 /iJy using different selection criteria. Upper 
limits are assumed considering S 24,1m = 200 fily. The distribution for the 
^-corrected data using a M82-like template is plotted by a thick black line. 



sources without photometric redshifts (compare thick continuum 
and dashed lines in Figure[5j does not introduce significant changes 
in the mean and scatter of the q 24 distribution. 

The fc-corrected mean and scatter values (biweight estimator) 
for the distribution of sources with S1.4GHZ > 300yuJy, mid-IR 
detections (> 200 /jJy) and estimated redshifts, changes to g 2 4 = 
0.71 + 0.47, 1.25 ± 0.47, and 0.77 ± 0.48 after using the M82- 
like, Arp 220-like and Lagache et al. model, respectively. The bot- 
tom panel in Figure [2] shows the M82-like fc-corrected ^24 values 
as a function of redshift for this sub-sample. It suggests a slight 
descending trend, from z < 0.5 to higher redshifts, probably asso- 
ciated with the presence of steeper mid-IR spectra for brighter ra- 
dio sources, although it is not possible to statistically discriminate 
a better representative SED between the tested templates. 

A simple linear regression to the M82-like /V-corrected g 2 4 dis- 
tribution, based on radio sources brighter than 300yidy and detected 
at 24pm, is given by q 2 dz) = (0.71 ± 0.04) + (-0.03 + 0.03) z. 
For this complete radio sample, we also find a flat distribution at 
lcr level, suggesting the validity of the FIR/radio correlation up to 
z ~ 3.5. Note, however that even at these bright radio flux densi- 
ties, at z > 1 the upper limits may overlap with the mean value 
found for the correlation, suggesting the possibility of confusion 
at higher redshifts between radio-loud AGN and those affected by 
mid-IR features in the SEDs. Future studies with submm and FIR 
facilities, such as SCUBA-2 and the Herschel Space Observatory, 



will help to better estimate the total bolometric FIR luminosity, and 
therefore discriminate between those sources affected by SED fea- 
tures or due to powerful radio AGN activity. 



4 RADIO-LOUD AGN ACTIVITY 

4.1 The transition form radio-loud AGN to star-forming 
dominated galaxies 

Making use of the M82-like fc-corrected data presented in the last 
section, we have worked out a method of selecting radio-loud AGN 
from our radio sample, based on different q 2 \ thresholds. 

It is well known that most powerful radio sources (> 1 mjy) 
are radio-loud AGN (Condon lll992T) . but the number counts ob- 
tained from deep radio surveys have suggested the appearance of a 
new, dominant population at fainter radio flux densities. It is well 
established that this faint radio population may be composed of 
star-forming galaxies (Condo rtl 19841) and radio-quiet AGN (Jarvis 
& Rawlings 2004). However, the transition from radio-loud AGN 
to star-forming galaxies is still not well defined. 

To look at the relative fractions of these populations as a 
function of radio flux density, we use the previous M82-like k- 
corrected data to identify radio-loud AGN based on different ^24 
cut-off thresholds (Donley et al. 2005). The main assumptions for 
the following results are the mean and scatter of the 24 yum/ 1 .4 GHz 
correlation (§3.3), and the flux thresholds from the radio and mid- 
IR maps. 

Figure[6]shows three different diagrams for the expected frac- 
tion of starburst/radio-quiet AGN ('star-forming-dominated' sys- 
tems0) and radio-loud AGN as a function of the 1.4-GHz flux den- 
sity. A value of q 2 4 = {q 2 4) - 2<x 924 (i.e. 2 a below the mean value) 
is adopted as the cut-off value of ^24 used to separate and dis- 
tinguish 'star-forming dominated' sources from radio-loud AGN. 
In particular, the cut-off used in the main panel of Figure [6] is at 
<?24 = -0.23, which in terms of luminosity means that radio-loud 
AGN are selected by having a radio power ~ 9 times or larger 
than that predicted by the mean value. All radio sources with mea- 
sured 24-ytim flux densities and ^24 values above the threshold are 



2 The reason we call these objects as star-forming-dominated sources 
comes only from the fact that they follow the FIR/radio correlation. Al- 
though, this is not entirely true because the radio emission from radio-quiet 
AGN (luminous Seyferts; sources that have been found to follow the cor- 
relation) takes the form of jets and suffers the same orientation-dependent 
beaming effects as the radio-loud AGN, implying it is most definitely AGN- 
related synchrotron emission. 
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Figure 6. The sub-mJy radio flux transition from radio-excess to starburst- 
dominated sources based on different 1724 thresholds using the M82-like 
^-corrected data. Sources without redshift have been included using a mean 
M82-like ^-correction (+0.33), non varying considerable the results but 
improving error bars. The "Radio-excess" fraction is composed of those 
sources detected at 24 fim, or with upper limits, with 1724 less than 2 <r from 
the mean value of the correlation. The "Star-forming Dominated" fraction 
is, therefore, that with ^24 values above the 2 <x threshold for radio-excess 
selection. The "Probable Star-forming Dominated" fraction is based on an 
estimation which corrects the incompleteness at 24//m of the "Star-forming 
Dominated" sources (see text for details). Hence, all the rest of the sources 
are very unlikely to follow the 24 fim/ 1 .4 GHz correlation and are called the 
"Probable Radio-excess" fraction. Error bars, based on the number counts, 
are 1-tr confidence limits given by Gehrels 1 19865). The main upper panel 
shows our best estimate of the different populations, based on the results 
from §3.3. The small panels at the bottom consider the same grey-colour 
criteria. On the left-hand side we modify the 524 mean and scatter based on 
the higher end of the correlation ((#24) + QVm ~ 1-2), and on the right-hand 
side we use the previous Appleton et al. 2004) results. 



We can then use this to correct the observed fraction of star-forming 
sources for incompleteness. In order to do so, we ran Monte-Carlo 
simulations for each radio flux density bin (from Figure [6]l, using 
the mean and scatter of the 24yum/1.4GHz correlation, to obtain 
the probability of detection based on the 24-fim image threshold 
at 200 /Jy. The correction for incompleteness to the "Star-forming 
Dominated" population is called, in Figure(6] the "Estimated Star- 
forming Dominated" fraction, which is mainly introduced at faint 
radio fluxes. The reminder of the uncertain population is classified 
as the "Probable Radio-excess" fraction. 

The main panel of Figure [6] shows our best prediction for 
the transition from radio-loud-AGN-dominated to a star-forming- 
dominated population at sub-mJy radio fluxes. It makes use of 
the mean and scatter of the M82-like fc-corrected ^24 values (from 
§3.3) to define the threshold value and the inputs to the Monte 
Carlo simulations. It is found that for star-forming sources with 
S i.4GHz = 200 - 500 /Jy (10~ 3 - 7 - 10~ 33 Jy), or with higher radio 
fluxes, the expected fraction of counterparts at 24 fim starts to be- 
come complete. Actually, the Monte-Carlo estimation predicts that 
> 85 per cent of the sources should have counterparts at 24 fim 
in this b in (it changes to > 98 per cent considering the Apple- 
ton et al j|2004l ^24 values). This confirms the previous estimation 
for the complete star-forming radio sample based on a cut-off at 
S i.4ghz = 300 //Jy. 

In order to show how the results change with uncertainties in 
the FIR/radio correlation, the two small panels at the bottom in 
Figure [6] use a n exag gerated scattering for the correlation, and the 
Appleton et al. 2004) results (left and right hand side, respectively), 
as indicated by the values above the panels. Each case, however, 
keeps almost the same fit to the high end of the q 2 4 distribution 
«<724> + cr 924 as 1.2). 

From our best fit we predict a significant (25 ± 5 per cent) 
fraction of radio-loud AGN at flux levels of 5 i.4ghz ~ 300 /Jy. The 
number of radio-loud AGN clearly increases as a function of the ra- 
dio flux density and above S i.4ghz ^ 3.2 mJy (10~ 25 Jy) essentially 
all radio sources are radio-loud AGN. At ~ 50/Jy flux density lev- 
els, the fraction of radio-loud AGN significantly decreases to zero 
(< 20% from the lower panel estimations), leaving the star-forming 
dominated systems as the only population at this faint radio flux 
regime. 

These simple analyses suggest that a large number of radio- 
loud AGN could be selected using deep mid-IR and radio sur- 
veys based on a simple radio-excess selection criteria. In terms of 
the AGN population, we note that our predictions for the "Radio- 
excess" fraction is just part of the total number of AGN, be- 
cause many radio-quiet AGN may be included in the 'star-forming- 
dominated' population. 



indicated in Figure [6] as 'Star-forming Dominated'. The so-called 
'Radio-excess' fraction corresponds to all those sources detected at 
24 fim with q 2 4 below the cut-off value, or undetected but with up- 
per limits below the threshold. The remaining population of sources 
will be those that are undetected at 24//m but have upper limits in 
g 2 4 that lie above the g 2 4 threshold for AGN selection. Those will 
be a mixture of starburst/radio-quiet systems and radio-loud AGN. 

We have already said in Section [3] that the faint radio popula- 
tion (< 300 /Jy) is affected by incompleteness at 24 fim. In fact, for 
sources that follow the FIR/radio correlation then, using the thresh- 
old for detection in the SWIRE data, we can predict the expected 
fraction of detections at 24/im as a function of radio flux density. 



4.2 Host galaxy properties 

To explore the optical galaxy host properties of the different popu- 
lations presented in §4.1, we estimated the rest-frame (U-B) colour 
(limited to sources with estimated redshift) for each source. The 
colour histograms for the three populations, "Star-forming Domi- 
nated" sources, the "Radio-excess" sources, and "All others" ("Es- 
timated Radio-excess" and "Estimated Star-forming Dominated" 
fractions combined), are plotted in Figure|7] 

Considering only the population with radio flux densities 
S 1.4GHz > 300 /Jy (left-hand panel), we find that both the "Radio- 
Excess" and "All others" populations present clearly redder dis- 
tributions compared to the the "Star-forming Dominated" popula- 
tion. A Kolmogorov-Smirnov (KS) rejects the hypotheses that ei- 
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Figure 7. Histogram of rest-frame U-B colours of all sources with es- 
timated redshift in two bins of radio density flux. Dashed, continuum, and 
dotted lines show the distribution of the "Star-forming Dominated", "Radio- 
excess" and "All others" ("Estimated Star-forming" and "Estimated Radio- 
excess") fractions, respectively. These correspond to the same ones plotted 
in the main panel of Figure|6] It shows that the Radio-Excess population has 
a fundamentally different nature to the Star-forming Dominated sources, 
and also confirms that at bright radio flux densities (> 300/jJy) the "All 
Others" population is comprised almost entirely of radio-excess objects. 



ther the "Radio-Excess" or "All others" sources are drawn from the 
same parent population as the "Star-forming Dominated" sources 
at > 99.999% significance in both cases, proving that they have a 
different nature. On the other hand, the colour distributions of the 
"Radio-Excess" and "All others" populations are statistically indis- 
tinguishable. This similarity is expected because, as argued in §3.3, 
at these radio flux densities the "All others" population ought to 
be composed essentially entirely of "Radio-Excess" sources (see 
Figure |6). At lower radio flux densities, the "All others" popula- 
tion should contain a mixture of starburst-dominated galaxies and 
radio-loud AGN, and the right-hand panel of Figure [7] shows that 
indeed at these lower radio flux densities this population has a more 
pronounced tail towards bluer U-B colours. 

The red colours found for the radio-excess sources may come 
from the absence of star formation or due to dust obscuration in 
these sources. Since they do not have detections at 24 pm - proba- 
ble absence of re-emission - we suggest that this radio-loud AGN 
population is hosted in red-galaxies with little or no star-formation. 
This is potentially a very useful way to detect obscured AGN 
missed in deep X-Ray surveys, which are res ponsib le of the bulk 



of the Cosmic X-Ray Background (Ueda et 



5 DISCUSSION 



esponstbl 
al jE003h . 



Previous studies of the universality of the FIR/radio correlation at 
high re dshift have combined 1.4-GHz data with ISO 15-pm data 
(Garr et ll2"002l) and Spitzer 24- and 70-pm data (Appleton et al. 
2004), both showing evidenc e of a correlation up to z ~ 1. More 
recently, Kovacs et al. 1120061) and Beelen et al. 1 2006) have shown 
evidence for the correlation at higher redshifts, for a sample of 15 
submm galaxies (z ~ 0.5 - 3.5) and for a handful of distant quasars 
(z ~ 2 - 6), respectively. In this work, despite the difficulty of mea- 
suring a reliable bolometric FIR emission, the monochromatic 24- 
pm flux densities show a tight correlation with the 1.4-GHz flux 
density up to z ~ 3.5. /c-corrections based on an M82-like template 



indicate that both the mean and the scatter of the ^24 values remain 
constant as a function of redshift (see Figure[3j. 

The mean and scatter of the distribution of q 2 4 values obtained 
after a M82-like /c-correction for all radio sources detected at 24 pm 
is <j>24 = 0.99 ± 0.32. This is in broad agreement with the results of 
Appleton et al. II2004I) . q 2 4 = 0.84 ± 0.28, based on a sample z < 1. 
However, if our sample is restricted to brighter radio sources for 
which all star-forming sources would be expected to have been de- 
tected at 24 pm, then the mean relation becomes q 2 4 = 0.71 + 0.47, 
with a scatter about twice that of Appleton et al. The median er- 
ror value for the ^24 's has been found to be about 0.1 dex, which 
suggests that this scatter is intrinsic for the correlation and does not 
come from observational errors. We can conceive of two possible 
explanations for this increased scatter at higher radio flux densi- 
ties. First, it may be that, because of the smaller number of sources, 
some radio-excess objects may fail to be excluded by the biweight 
estimator, and therefore both bias and increase the scatter of the 
relation. Second, it may be that for these more powerful sources, 
the M82-like template (a relatively low luminosity source) is not 
an accurate representation of the FIR emission of all sources, but 
that a broader range of SED types is present, increasing the scatter. 
Figure [2] provides some evidence for this: although a very flat q 2 4 
distribution is found in the full radio sample (top panel), some red- 
shift variation may be present in the brighter radio sources (lower 
panel), in particular the presence of a handful of sources with higher 
^24 values at low redshifts z < 0.5, which would be more consistent 
with a steeper mid-IR slope, such as that given by the Arp220-like 
template. 

Recently, Boyle et al. M2007I) have found a significantly higher 
mean value, q 24 = 1.39 ± 0.02 (2 a away f rom o ur estimation) 
whilst applying stacking techniques (Wals et al 2005) to ATC A 1.4- 
GHz images using positions from Spitzer detections in the CDFS 
and ELAIS fields. Their work lacks redshift information, however, 
and therefore does not include /c-corrections. We have shown the 
relevance of different mid-IR features for high redshift sources 
at the flux regimes of this work, suggesting a clear underestima- 
tion for the intrinsic 24pm/1.4GHz scatter in their work. Also, a 
stacking analysis for a mixture of star-forming galaxies and radio- 
loud AGN does not discriminate an exclusive correlation for star- 
forming galaxies. Since the flux thresholds from our work are sim- 
ilar to their ones, we suggest that selection effects based on differ- 
ent sample selection criteria (from radio or mid-IR detections) may 
also be affecting the scatter on the 24/jm/1.4GHz correlation. 

The FIR/radio correlation in not precisely linear, specially in 
optically selected samples containing faint galaxies. It is expected 
that the radio emission, in low luminosity galaxies, is more defi- 
cient since cosmic rays a re more likely to escape from the galaxy 
(Chi & Worfendal Jl990T) . At redshifts larger than unity the Cosmic 
Micro wave Background might quench the radio emission (Condon 
II 1992b and a cosmic evolution for the magnetic fields cannot be 
rejected too. Nevertheless, in this work the exploration of more de- 
tailed models, showing for example luminosity (or mass) depen- 
dencies for the 1724 parameter, was limited by the uncertainties in 
photometric redshifts. 

A simple radio flux estimation (S'"Iq Hz as 35//Iy) shows that 
a starburst galaxy with a typical M82 luminosity can be observed 
up to moderate redshifts only, z ~ 0.4. In fact, it is expected that 
the observed high-redshift high-luminosity Universe is composed 
mainly of ULIRGs and AGN. U LIRGs are usually related to in- 
teracting systems (Sanders et al. ll 19881) and are expected to suf- 
fer strong absorption, implying large ^-corrections for the 24-pm 
emission. Insights for steeper mid-IR SED at Si.4ghz > 300yuly 
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are hinted (bottom Figure [2j, however, the lack of a clear scatter 
excess in the FIR/radio distribution at z ~ 1-2 (Figure Q} sug- 
gests that the bulk of the ULIRG-like high-redshift sources from 
our sample do not suffer large variations due to silicate absorp- 
tion. The correlation found in this work is therefore quite unex- 
pected, but suggests that a M82-like mid-IR template, intrinsically 
coming from a lower luminosity source, is still a good representa- 
tion for the mid-IR SED of high-redshift star-forming-dominated 
systems. The M82-like template is also supported by Lagache's 
ULIRG (L IR = 10 I2 M G ) model, and particularly leaves Arp220 
as single extreme c ase. Th ese results are more relevant considering 
that Kasliwal et al. H2005I) estimate that about half of ULIRGs at 
z ~ 1-1.8 may be unobserved at 24 fjm due to silicate absorption at 
9.7 yum in the Bootes field, a suggestion not particularly supported 
by our data. For ULIRGs selected in sub-mm observations, their 
mid-IR SED is still under debate, bu t temp lates such as M 82 and 
Arp 220 (Menendez-Delmestre et al. ll2007l and Pope et al ]|2006l. 
respectively), have been found as good fits for their mid-IR range. 

In Section [4] we have shown the capabilities of deep mid-IR 
and radio imaging for the selection of radio-loud AGN ("radio- 
excess") based on different cut-offs in q 2 4- We have shown that at 
S 1.4GHz > 300 yuJy, selected radio-excess sources strongly suggests 
nuclear activity. The radio-excess sources present a redder distribu- 
tion of rest-frame U-B colours compared with those identified as 
star-forming-dominated galaxies. This strongly suggest a different 
nature between the populations selected by the simple q 24 threshold 
at -0. 23. Previous population synthesis models (Jarvis & Rawlings 
2004) have predicted that radio-excess s ources at these flux den- 
sities are mainly FRI (Fanaroff & Rilev ll 19741) . It is known that 
FRI d o not f all into the standard unified scheme for AGN (An- 
tonucci HT993h and show no bright AGN nucleus. Therefore their 
red colours indicate that these radio-loud sources do not present 
significant star-formation activity and/or are hosted in red massive 
galaxies. 

Jarvis & Rawlings M2004I) predict that approximately half of 
the sources with radio fluxes between lOOyuJy =S Si4GHz < 300yuJy 
present radio-loud activity. We have predicted from our best fit a 
fraction of 25 ± 5 % at this flux regime, in disagreement by a factor 
of two with their estimation. Figure|6]shows that Appleton et al. es- 
timation may explain a 50% fraction, although we demonstrated in 
§3.3 that incompleteness biases the assump tions of those results. 
On the other hand, Simpson et al. |j|2006) tentatively suggested 
a radio-quiet AGN fraction of 20% at this radio fluxes regime, 
but since radio-quiet sources have been found to also follow the 
FIR/radio correlation, in this work we are not able to discrimi- 
nate between this population and starburst galaxies. Therefore, this 
radio-quiet AGN population remains uncertain. 



6 CONCLUSIONS 

We have analysed the 24 /im properties of a radio-selected sample 
in order to explore the well-known FIR/radio correlation at high 
redshifts. We have shown evidence that 24-/jm data do not give 
the best estimation of the total bolometric FIR emission, reflected 
in a considerably larger scatter on the correlation than p revious 
work using longer wavelengths (IRAS/60 //m - Yun et al. II2001L 
SpitzerPO (im - Appleton et al. ll2004l) . Nevertheless, despite the 
large scatter, the IR-24 yum/radio- 1 .4 GHz correlation seems to ex- 
tend up to redshift ~ 3.5 with roughly the same mean and scatter. 
The fact that we do not observe considerable variations in the cor- 
relation as a function of redshift, mainly seen on the well defined 



upper end of the /c-corrected q 2 4 distribution (see top Figure [2] & 
[3}, highly suggests its validity all the way back to primeval times. 
This is the first evidence for the extension of the correlation over 
this very wide range of redshifts. 

Statistically, the correlation is described by the basic observ- 
able parameter q 24 , the ratio between the flux densities at 24yum 
and 1 .4 GHz, respectively. The q 2 4 factor parameterise the link (in 
the local Universe at least) that star-forming galaxies have between 
the synchrotron emission, from cosmic rays accelerated by massive 
supernovae explosions, i.e. the massive end of the mass function in 
galaxies, with the reprocessed UV photons re-emitted by the dust 
at FIR wavelengths. Besides an intrinsic scatter for the correlation, 
the presence of an AGN may also modify the mid-IR spectra and 
the radio emission, biasing, and therefore increasing the scatter on 
a (j> 2 4-b ase d correlation. We found a mean value (biweight estima- 
tor) of q 24 = 0.30 ± 0.56 for the observed and q 24 = 0.71 ± 0.47 
for the fc-corrected (M82-like) data based on an expected complete 
radio sample (S 1 4ghz > 300 yuJy) at 24 ^m. 

Three different /^-corrections have been tested: two based on 
well-known local starburst sources (M82 and Arp 220) and one on 
a SED ULIRG model (Lagache et al. ll2004l) . The main difference 
introduced by using different templates is caused by the silicate 
absorption feature at ~10yum being redshifted into the observed 24- 
fim band at z ~ 1-2 (Figure [TJ. Based on the smallest scatter 
found for q u data, we propose an acceptable /t-correction given by 
an M82-like mid-IR template, which suggests no strong evidence 
for extreme silicate absorption in the bulk of our sample. We do, 
however, find tentative evidence for a steeper mid-IR spectrum than 
M82 for brighter radio sources > 300/^Jy, in order to explain the 
slight offset from the correlation for sources with z < 0.5. 

We have found a clear dependency for the detected num- 
ber of star-forming-dominated galaxies (starburst and radio-quiet 
AGN) as a function of the radio flux density. Making use of the 
derived 24yum/1.4GHz correlation, and the flux density threshold 
from the 24-yum image, we ran Monte-Carlo simulations to esti- 
mate the incompleteness at 24yum, and therefore the intrinsic frac- 
tion of star-forming-dominated galaxies. This allowed us to esti- 
mate the probable fraction of radio-loud AGN and to describe the 
transition from AGN-dominated to star-forming-dominated galax- 
ies at faint radio flux densities (< 1 mjy). Our best-fit estimation 
predicts a significant fraction (25 ± 5 per cent) of radio-loud AGN 
at S 1.4GHz ~ 300 yuJy, rising quickly to become a dominant AGN 
population at brighter radio flux densities (> 1 mJy), and falling to- 
wards zero (certainly < 20 per cent) at fainter fluxes. The se resu lts 
are in broad agreement with previous studies (e.g. Condorjl984). 

We have also shown that the rest-frame U-B colours, of the 
sources that are expected to be radio-loud AGN, have a redder dis- 
tribution than the starburst/radio-quiet AGN galaxies (Figure [7J, 
strongly suggesting a different nature between the sources selected 
by our simple q 24 = -0.23 threshold. These radio-excess sources 
seem to be hosted in red galaxies with little or no star-formation im- 
plied by the non-detection at 24yum. This is therefore an extremely 
prom ising way to select obscured Type-2 AGN (Comastri et al. 
1 19931) . with a radio-loud nature, which are missed by deep X-ray 
observations. Nevertheless, future spectroscopic follow-up of these 
sources, and also a better description of the FIR emission, is needed 
to fully test this method. 
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